EVOLUTION, INTERACTION AND COLLISIONS 
OF VORTEX FILAMENTS 
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Abstract. Several progresses have been done very recently on models for the 
dynamics of one or more vortex filaments in 3-D fluids. In this article we survey 
the recent and previous results in this topic. We finally present a collection of 
new situations of filaments collapse. 
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1. Introduction 



The purpose of this work is to describe some aspects of the dynamics of ho- 
mogeneous three-dimensional incompressible fluids. The evolution of such fluids is 
governed by the Euler equations 

^ dtuj + V ■ Wuj = uj ■ Wv, divw = 0, (1.1) 

where : K x — >• and u — curl(w) : M x M'^ — > denote the velocity and the 
¥^ vorticity of the fluid. 

k>< We shall focus on fluids in which remarkable structures, called vortex fllaments, 

5_] are present. A vortex fllament is a very thin tube in which the vorticity is sharply 

concentrated. It can be asymptotically assimilated to a curve s) in M'^, where 
< € M denotes the time and s G M the arc-length parameter. According to a suit- 
able formal derivation, which will be described in the next section, the asymptotic 
motion law of the curve is the binormal flow 

Xt=dsXAd^,X- (BF) 



Equation (BF I exhibits a rich variety of motions, which will be presented in Section 
[2] below. Particular attention will be payed to motions generating a singularity in 
finite time. 

On the other hand, special solutions of the binormal flow are given by the infi- 
nite straight filaments. Up to a rotation, we can assume that there are parallel to 
the ea-axis so that they can be parametrized as xC^j"') — i^Oi<^) for some fixed 
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Xq £ E^, where cr G M can be also seen as an arc- length parameter. Of interest is 
the evolution of perturbations of such straight filaments, namely the nearly paral- 
lel vortex filaments, which are parametrized by x(t, cr) = (Xq + ^(t,(T),(T) where 



^'(t,cr) S M ~ C is small. It turns out that linearization about (BF| leads to the 
linear one-dimensional Schrodinger equation for 

iat* + a^* = 0. (1.2) 

Next, we wish to investigate the dynamics of several vortex filaments that are all 
nearly parallel to the same ea-axis. We consider a collection of iV > 2 filaments, 
which are represented by the complex- valued functions (t) € C, with i € K 

and (T e M. In 1995, Klein, Majda and Damodaran [39] derived, in a suitable 
asymptotic regime, a system of coupled equations for the filament positions ^/j 

idt^, + a,r,dl^, + ^ Ffe - 0, 1<J<N. (1.3) 

Here Tj € K corresponds to the circulation of the j-th filament, and the parameter 
aj € M is related to the core structure of the filament. In particular, in the case 
where the vortex filaments are all exactly parallel, i.e. if *I'j(t, cr) — Xj{t) does not 
depend on a, the previous system reduces to a two-dimensional system of ordinary 
differential equations for the positions Xj{t) 

dX,{ V — A _ X^ — Xu , , 

This system is called point vortex system or Kirchhoff law, and it has been inten- 
sively studied in the literature. 



We observe that the system (1.3) combines two different aspects of the dynam- 



ics. On the one hand, the linearized self-induced motion of each filament -already 



expressed in the equation (1.2)- is represented through the Schrodinger operator. 
On the other hand, the interaction of the j-th filament with the other filaments is 
represented by the potential field X^fc^^j !"*.•(* j ~^k)l\^ j — '^k^ ^ which corresponds 
to the velocity generated by the other filamentf[^ For straight filaments there is no 
self-induced motion, which means that each filament moves only with the velocity 
induced by the other filaments according to ( |1.4[ ) - exactly as for planar vortices, 
see Section |3] hereafter. The fact that both effects are taken into account in (1.3) is 



really due to the asymptotic conditions chosen for the derivation, which relate the 
wavelength of perturbations, the core sizes and the separation distances between 
the vortex tubes in a particular way. 



Since the system (1.3) is thought to be a simplified model for the dynamics of 
filaments, a first natural issue is the existence and uniqueness of solutions. The 
potential fields are not well-defined at points when two or more filaments collide 
in finite time. One possibility is to consider weak distributional solutions, allowing 
possibly for exceptional collisions, and for which the potential is well-defined almost- 
everywhere (see Section |4|. But then uniqueness is not achieved in this class. 
Therefore, we will focus on a class of stronger solutions, for which the filaments 
are well-separated up to the largest time of existence of the system. Hence the 
largest time of existence corresponds to the first collision time in this framework. 
It should be mentioned that, even for the simpler situation of straight filaments (the 
point vortex system), collisions in finite time are known to occur for some initial 
configurations (see Section [3| . Therefore we cannot hope for global existence in 
general. We shall determine sufficient conditions on the initial collection of vortex 



^According to the Biot-Savart law, see Section jS] 
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filaments that lead to large time or global existence, and we shall also study the 
case of finite-time collapses. 

The remainder of this work is organized as follows. In Section [2] we present the 



formal derivation of the binormal flow ( BF ) for the motion of one single filament and 



we describe some of its remarkable properties. In particular we review the recent 



results on singularities for self-similar and almost self-similar solutions of (BF), as 



well as for their perturbations. Section|3]contains a brief state-of-the-art concerning 



the point vortex system (1.4), as an introduction to the more complicated dynamics 
of almost parallel vortex filaments. Then, Section [4] gives an overview of existence 
and uniqueness results concerning the system ( |1.3| . In particular, we present some 
recent results from |6j on symmetric configurations of vortex filaments. We conclude 
by providing several and simple examples of finite-time collapses of vortex filaments. 
The collapses are of different nature, and some of them are new. 

Notations. In the following we will systematically identify complex numbers and 
real vectors. In particular, introducing the matrix J = we shall identify 

the M^-vcctor JX and the complex iX. 

Acknowledgements. The first author is partially supported by the ANR project 
"R.A.S." . Both authors are grateful to Vincent Torri for his support in using Scilab. 

2. One filament 

In this section we briefly present classical facts about the local induction ap- 
proximation and about the binormal flow, as well as recent results of singularity 
formation in finite time. For parts of this description the reader may also consult 
the book 52J and the surveys [62 and JTP. 

2.1. The Local Induction Approximation. We shall start with a short descrip- 
tion of the way the binomial flow equation was flnd as a model for vortex filaments 
dynamics in a three-dimensional ideal incompressible fluid. A fluid with a vortex 
filament is a fluid where the vorticity uj(t) is a singular measure supported along 
a curve which is parametrized by arc-length parameter s, and with density 

rdsx{t)- Here F denotes the constant circulation along the fllamcnt. As a conse- 
quence of Kelvin's circulation law vortex filaments move with the flow, so in order 
to get the dynamics of the vortex curve we have to compute the velocity of 
the fluid near the curve. The three-dimensional Biot-Savart kernel is . ] , so the 
velocity of the fluid is 

/ 1 \ r f°° X - y(t s) 

Vit, = -V A * .(t, y_ |.-;,(M)|3 ^ 

We fix t, we suppose without loss of generality that x(*)0) = (0) 0) 0), 9sx(^, 0) = 
(0, 0, 1) and we get focused on what is happening for the velocity of the fluid near 
(0,0,0). For this purpose, two localizations are done, one is that x(^: •s) will be 
approximated by a Taylor development of order 2 near s = and the other is that 
the integral has to be considered only locally around s = 0, on a segment [— L, L]. 
We obtain that v{t, (e, 0, 0)) is approximated by 

F (e,0,0)A(0,0,l) ds ^ T dsx{t,0)Adh{t,0) 



AtT £2 7_^|l + s2|5 AtT 2 7-^|l + s2|S 

When e -> 0, the first term is diverging but in the same way the velocity does 
around one straight vortex filament. One straight vortex filament remains still in 
a fluid, so this first contribution should not influence the dynamics of one single 
vortex fllament Finally, even after localization, the integral in the last term is 
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still diverging as | lne|. A time rescaling allow to resorb this growth and then the 
binormal flow 

dtX^dsXAd^X (2.1) 
appears as an approximation for the dynamics of vortex filaments, called Local In- 
duction Approximation. It was discovered by Da Rios 21J in 1906 and re-lightened 



by Arms and Hama [5] in 1965. However the rigorous derivation of (2.1) from the 
Euler equations ( |1.1[ ), and in fact even the question of the stability of the vortex 
filaments, are open. 

Being a crude model, the binormal flow has advantages and disadvantages. The 
positive points are the fact that special solutions as straight lines, translating circles, 
helices are on the one hand solutions of the binormal flow and on the other hand 
these dynamics are observed in fluid dynamics. This way other special solutions of 
the binormal flow may give some intuition for particular type of vortex filaments 
in fluids. A successful example in this direction are the traveling type solutions of 
the binormal flow found in |32j and displayed in an fluid experiment |33) . Also, 
the binormal flow is accessible for numerics. One of the main issues of the L.I. A. is 
that closed curves conserve their lengths under the binormal flow which is not the 
case for closed vortex filaments in fluids. 

The binormal curvature flow also appears in the context of supcrfluids that are 
governed by the Gross-Pitaevskii equation 

idt'^ + A^ + ^{l-\^\^) = 0. (GP) 



More precisely, (BF ) is conjectured to govern the asymptotic dynamics of codimension- 



2 submanifolds around which the Jacobians of the solutions to (a rescaled version 



of) ( GP ) concentrate. This was proved in any dimension N in the special case of the 
N — 2-dimensional sphere [53] . In a parallel setting, mean curvature flow governs 
the vortex dynamics for parabolic Ginzburg-Landau equation [T3]. The fact that 
vortex filaments exhibit a common asymptotic motion law in fluids and superfluids 



results from a fundamental analogy between ( GP ) and the Euler equations in such 
regimes. More precisely, the Madelung [55] transform ^' = y/pexp{i(p), as long 
as ^ does not vanish, yields a compressible Euler type system for the variables p 



and V = S7ip. It is called hydrodynamical form of (GP) and it reduces formally 



to the incompressible Euler equations for v in the above-mentioned singular limit 



of (GP). Further details may be found in the recent survey and in references 



quoted therein. 

2.2. The binormal flow. The binormal flow is a rich geometric equation. It is a 
completely integrable system. In the following we shall denote curvature and torsion 



of a solution of (2.1) by c{t,s) and T{t,s). Some remarkable conserved quantities 
of the binormal flow are the kinetic energy J c^(t, s)ds, helicity / c^T{t, s)ds, linear 
momentum J X ^ dsXi^j s)ds, angular momentum / X A (% A dsx)it, s)ds. It is 
a reversible in time equation, invariant with respect to translations and rotations. 
The tangent vector T{t, s) of a solution of the binormal flow satisfies the Schrodinger 
map equation 

dtT^TA dlT (2.2) 

which plays an important role in ferromagnetism as a simplification of the Landau- 
Lifschitz equation. 

Let us present now the link with the Schrodinger equation. As mentioned in 
the introduction, the linear Schrodinger equation can be obtained easily as a rough 



approximation of the binormal flow. This can be seen by plugging in (2.1) the 
ansatz 

X(i,s) = {ex{t,s),ey{t,s),s). 
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To leading order when e — >■ 0, we obtain that x + iy is a solution of the linear 
Schrodinger equation. It was in 1972 that Hasimoto |32] has discovered the link 
between the binormal flow and the cubic one-dimensional nonlinear Schrodinger 
equation. The so-called filament function constructed from curvature and torsion 



of a solution of (2.1 1 



satisfies the 1-D cubic NLS 

idt^ + + \ (IV-P - A{t)) ^ = 0, 

with A{t) in terms of curvature and torsion (c, T)(t, 0). More precisely, 

, dsc^ -ct"^ . 2 



A(t)= (^2^^^— +c^j (i,0), 

which might be seen as a compatibility condition. Let us notice that the Hasimoto 
transform can be seen as an inverse of Madelung's transform - actually curvature 
and torsion satisfy an Euler with quantum pressure type of system. One can remark 
that the non- vanishing of the curvature appears as a constraint. Nevertheless, Koiso 
gave in [iD] a way to avoid this issue by using, instead of Frenet frames, parallel 
frames (T, 61,62) that satisfy 

( T \ / a p \ ( T 

ei = -a ei 

V 62 / V -/3 / V 62 

Then a -I- z/3 is a solution of 1-D cubic NLS. For details of this kind of Hasimoto 
transform we refer to Appendix A of 62 . As a conclusion, up to a change of phase, 
the filament function satisfies the 1-D cubic NLS. 

By using the Hasimoto transform the first local well-posedness results were ob- 
tained for the binormal flow - curvature and torsion were considered in high order 
Sobolev spaces, see BOIES] . 

A recent result of local well-posedness for the binormal flow, for less regular 
closed curves, was obtained by Jerrard and Smets [35|, I36j by considering a weak 
version of the binormal flow. Their method acts at the level of the tangent of the 
curve and it does not use the Hasimoto transform. Also, they have derived this 



way new results for the Schrodinger map equation (2.2). 

As far as we have seen, there exists a number of intimate relations between the 
evolution of filaments in fluids and dispersive equations. We end this subsection 
by discussing other similar connections. In order to improve the validity of the 
binormal flow as a model for vortex filaments dynamics, an extended version of 



(2.1) was considered by Fukumoto and Miyazaki [23], which includes an axial flow 



along the filament, 

dtX = dsX A a,'x + a ( d'tx + ^d'^x A {dsX A d'ix) 



Here a denotes the magnitude of the axial flow. The Hasimoto transform leads 
then to the so-called Hirota equation 

which is in link both with NLS and the modified KdV equation (see [46] for the 
initial derivation of this equation). It actually turned out that the modified KdV 
equation appears as a model for the evolution of the curvature of a vortex patch in 
a two-dimensional fluid [551 ES] , see also [30] for more details and for self-similar 
solutions of this model. 
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2.3. Singularity formation scenarios for the binormal flow. An interesting 
class of solutions that is not covered by the previous well-posedness results are the 
self-similar solutions of the binormal flow and their small and smooth perturbations. 



Since solutions of (2.1 1 are considered to be parametrized by arc- length, self-similar 



solutions are of type 



It is then easy to compute that they form a family of solutions Xa indexed by a 
real parameter a > and {ca,Ta){t, s) = (^-^, ^) (BS])- They have been studied 
in the last decades as a frictionless model for the evolution after a line-line recon- 
nection in the absence of counterflow for superfluid 4He and in ferromagnetism 
[531 [ini Sni SB SI] • They also appear, this time for finite length filaments with ap- 
propriate boundary conditions, in some models for the fiber architecture of aortic 
valve leaflets, see [61]. Their rigorous mathematical study was achieved by Gutier- 
rez, Rivas and Vega, who gave in |29j a description of the profile and in particular 
established that 

Xa(0,s) = sA+I[o,oo)(s) - S-4"I(_oo,0](s), 

with e distinct, non-opposite and sin ^^^'^^ ~ e^^. So self-similar 
solutions of the binormal flow are smooth infinite curves that generate a corner in 
finite time. In fluid mechanics, their dynamics is close to the one of the delta-wing 
vortex. In ^22^ numerical simulations for self-similar solutions of the binormal flow 
were shown to be in correlation with the physical experiment. 



Another family of (2.1) generating a singularity in finite time at the level of the 



Frenet frame are the almost self-similar solutions 



x{t,s)^e^'°^'V~tG(^-^^ 



where is a real 3x3 antisymmetric matrix, see [151 [50] . Gutierrez and Vega [3U] 
proved that 

Xa(0,s) = se^i°«I^U+Voo)(s) - se^>°sl^lA-I(_oo,o](s). 

The curves at times t — are spirals, some of them with a singularity point, see 
[68j for some numerical simulations. Both similar and almost self-similar solutions 
have self-similar curvature and torsion, so the filament function is also of type 

ip{t, s) — -^G and has to solve 

i5,^ + a2^+i(|V,|2_^)^^0. (2.3) 
Therefore the data has to be homogeneous of degree —1. The case of Dirac distribu- 



tion leads to the self-similar solutions of (2.1 1, and the principal value distribution 



to the almost self-similar solutions respectively. 

The study of small and smooth perturbations of self-similar solutions was initi- 
ated by Banica and Vega [TllHlinilin] from the point of view of the stability of the for- 
mation of the singularity. The starting point of the analysis is the Hasimoto trans- 
form: the filament function of the self-similar solutions Xa is ipa{t,x) = ^e**^. 



solution of (2.3 1. By the pseudo-conformal transformation 



^{t^x) = T{u + a){t, x) ^ — e'^Tt u + ai -,- 

V i \t t 
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understanding the behavior at time for perturbations of tpa, solutions of (2.3), is 
equivalent to understanding the large time behavior for the equation 

idtu + dlu + ^ (k + ap - a^) (w + a) = 0. (2.4) 

A scattering theory has been developed in fS'.'S] for this problem, by strongly relying 



on Fourier analysis. The informations on (2.4) have been transported at the level 
of the binormal flow by treating differently the regions ts^ < 1 and 1 < ts^ . In 
order to discard the linear non-oscillating term, the asymptotic profile for u{t) has 
to be modified into the long-range one 



; t itd: 



In [5j it has been shown that wave operators exist, with a strong H'' decay, 
for asymptotic profiles u+ S n H'^ n VF*'^. By adding the extra-conditions 



x^Uj^ £ I? and s > 3 the existence of solutions of (2.1 1 that develop a corner- type 
singularity at i = has been derived. This insures that the formation of singularity 
in finite time appearing in the self-similar case is not an isolated phenomenon. A 



similar result for the almost self-similar solutions of (2.1 1 has been proved in |31) . 

Asymptotic completeness with loss of regularity was then proved in [5| , in spaces 
defined in terms of the behavior of the low Fourier modes. At the level of the 
binormal flow, this information has been translated into the fact that all small and 
regular perturbations at time tg 7^ of a self-similar solution Xa still yield, after 
evolving by the binormal flow, a singularity at time t — Q. 

Finally, in flU] we give a geometrical description of the curves at the singularity 



time by getting a control of the evolution of weighted norms for equation (2.4). 
We prove the stability of the selfsimilar dynamics of small pertubations of a given 
selfsimilar solution, and in particular the fact that the angle of Xa is recovered at 
time t = 0. 



3. The point vortex system 



This section is devoted to a short presentation of the point vortex system (1.4). 
The reader may find the subsequent results as well as many additional details in 
the surveys [H [SHI H] and in references quoted therein. 

The point vortex system arises as an asymptotic motion law for point singu- 
larities in two-dimensional incompressible fiuids, which are governed by the incom- 
pressible Euler equations, or superfiuids, governed by the Gross-Pitaevskii equation 



(GP ). In contrast with the three-dimensional case, the persistence of point vortices 
and the validity of the point vortex model have been established rigorously for 
well-prepared data, at least up to the first collision time between the vortices (see 
[MIIMI for the Euler equations and, e.g., [20l[T3| for the Gross-Pitaevskii equation). 



System (1.4) possesses an Hamiltonian structure, since it has the form 



^A^A^) + JVx,i/(Xi, . . . ,XAr) = 0, 1 < J < iV, 



' dt 

with the Hamiltonian H defined by 



i/(Xi , . . . , Xjv) = - J ^ r.Tfe In \X, - Xk 



2 



As a consequence, we have the conservation of H 

H{Xi{t),..., XN{t)) = i/(Xi(0), . . . , XAr(O)). 
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Moreover, the point 

N N 

and the angular momentum 

N N 

are preserved in time. If Tj ^ 0, we define the constant center of vorticity 

AT 

c= (Er,)-ic. 

It follows that the quantity 

j,k j,k 

is constant as well. A direct consequence of the previous conservation laws is that, in 
case of circulations Tj all having the same sign, no collision between the vortices can 
occur in finite time and the A^-vortex system admits a unique and global solution. 

It is well-known that the point vortex system admits particular solutions that 
keep the geometric shape of the polygon formed by the vortex positions. To describe 
these solutions we shall adopt the following definitions from, e.g., [56]: we say that 
the configuration {Xj{t))j is 

• a fixed equilibrium if Xj{t) ~ Xj{0), Vj — 1, . . . , N; 

• rigidly translating if, for all t, there exists V £ C* such that Xj{t) = V, 
Vj = l,...,iV; 

• a relative equilibrium if there exists some Xq e C and for all t, there exists 
Lu eR* such that Xj{t) = iuj{t){Xj{t) - Xo), Vj = 1, . . . , A^; 

• a self-similar (or homographic) motion if there exists some Xq G C and for 
all t, there exists w S C, with 3m(a;) ^ 0, such that Xj{t) — iijj{Xj{t) — Xq), 
Vj = l,...,Ar. 

Note that in the first three cases both the shape and the size of the polygon 
formed by the points are preserved, since the motion is a translation in the second 
case and a rotation- with angular velocity uj{t)- in the second case. In the case 
of homographic motion the configuration at time t is obtained from the initial 
one by a rotation and a dilation. More precisely, it must satisfy Xj(t) — Xq = 
{Xj{0) — A'o)-\/T^^li7^cxp(iA In |1 — </io|), where the reals to and A are determined 
in terms of uj. When to > 0, which actually corresponds to the case 9m(w(t)) > 0, 
the configuration shrinks and collapses into the center of vorticity Xq at time to- 
Otherwise we have Iq < 0, i.e. 5m(a;(t)) < and the configuration expands in 
time. 

For N = 2 the distance |A"i(i) ~X2{t)\ remains constant therefore the solution is 
global. The motion is a uniform rotation about the center of vorticity if Fi -|-r2 ^ 0, 
and a rigid uniform translation otherwise. 

For A^ = 3 the motion is integrable and has been widely investigated, see e.g. 
[28l [H m |3] . An important part of the analysis was initiated in 1877 by Grobli 
[28] in his PhD thesis. In particular Grobli noticed that the 3- vortex system can 
be entirely formulated in terms of the three distances between the vortices. The 
motion can then be analyzed qualitatively through the use of trilinear coordinates 
and phase diagrams, see the work of Aref [T]. There exists a large variety of regimes: 
fixed equilibria (if Fir2 + Fir3 + r2r3 = 0), rigidly translating equilateral triangles 
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(if Fi + r2 + Fa =0 and T ^ 0), uniformly rotating equilateral triangles (if Fi + 
r2 ^0 and T 7^ 0), coUinear rotating or fixed configurations. Also, self-similar 
(shrinking or expanding) collapsing configurations exist pSl [1] [521 13j • It turns out 
that the conditions T = and Fir2 + TiT^ + T2T3 = are both necessary and 
sufficient to obtain a self-similar expanding or shrinking motiorj^ and this provides 
an effective and simple way to construct them explicitly. On the other hand, the 
linear stability of such motions has been investigated in, e.g., |38[ l3] . 

We shall come back to these collapses in Section |4.3[ where we will use them to 
construct collapsing filament solutions. 

For N > 4 the A^-vortcx problem is not integrable in general, but it still exhibits 
some remarkable solutions. First well-known relative equilibria are given by the 
uniformly rotating configuration of N identical vortices placed on the same line in 
a particular way, or, when N is odd, by — 1 identical vortices and a central vortex 
of arbitrary circulation. Other relative equilibria are obtained by setting one vortex 
with same circulation Tj = F at each vertex of a regular N - polygon. The motion 
is then a uniform rotation with constant angular velocity w = F(A^ — f)/(2i?^), 
where R denotes the size of the polygon. The problem of stability of these iV- 
polygons was raised by Thomson (Lord Kelvin) in the nineteenth century and has 
attracted much attention since then. The full answer to Thomson's question - that 
the regular A^-polygon is stable if and only if A < 7 - was finally recently provided 
by Kurakin and Yudovitch [42l |43] . One can get other uniformly rotating relative 
equilibria by adding to a regular A-polygon configuration one point of arbitrary 
circulation Fo at the center of the polygon, in which case the angular velocity is 
Lu = (F(A^- 1) + 2Fo)(2i?2). In particular, choosing Fq = -F(A^- l)/2 one obtains 
a fixed equilibrium solution. The stability of centered A^-polygons has been proved 
to hold by Cabral and Schmidt [T7] for all A' provided the circulation Fq lies in 
a suitable interval depending on A. In particular, the range of stability includes 
negative values of Fq for A^ < 7. 

Finally, we mention that other relative equilibria can be constructed in the form 
of concentric rings of polygons with identical vortices. The polygons may be not 
regular, and the rings may have different circulations (heterogeneous rings). Rel- 
ative equilibria of identical two-vortex or three- vortex rings have been determined 
by Aref, Newton, Stremler, Tokieda and Vainchtein [3], and O'Neil ;58J constructed 
relative equilibria of heterogeneous three- vortex rings. 

The existence of collapsing solutions for A^ > 4 is a difficult issue. Given a 
set of circulations such that J2 j^j for all J C {I, . . . , A^}, the set of initial 
configurations leading to a collapse in finite time has zero Lebesgue measure in M.^-'^ 
|53j . On the other hand, in general it is not known whether self-similar collapses 
exist. Necessary conditions for such motions are T = and FiF2-l-F2F3-|-F3ri = 0, 
but, in contrast with the case A^ = 3, there are not sufficient. Nevertheless, some 
specific cases have been explored: four and five- vortex collapses were studied by 
Novikov and Sedov [59] and collapses of two- vortex rings were constructed in [41]. 
Collapsing configurations for infinite point vortex lattices are constructed by O'Neil 
in [ST]. Finally, [SS] exhibits self-similar collapses of heterogeneous three- vortex 
rings. 



4. Several filaments 

4.1. The interaction of several filaments. The objective of this paragraph is 
to review known existence results concerning the system (1.31, which have been 
established by Klein, Majda and Damodaran [39] and Kenig, Ponce and Vega [37] . 



^If the initial configuration is not a fixed equilibrium. 
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We also mention the work by Zhakarov (SnUlZ]- We recall the system (1.3) supple- 
mented by initial data 5'j,o- 



0, 



1 < J < ^, 



4.1.1. Some formal properties. For any > 1 and any choice of circulations Tj G 
the dynamics formally preserves the following quantities: 
The Hamiltonian 

N 



j=i 



3=^k- 



The mean angular momentum 



/ 



N 



i=i 



and the mean center of vorticity 



N 

i=i 



They are the three-dimensional analogs of the notions of Hamiltonian, angular 
momentum and center of vorticity for point vortices, which have been defined in 
Section [Sj There is also conservation of the quantity 



r N 



which can be interpreted as a mean momentum. 

Finally, in the case where the vortex core parameters and the circulations satisfy 

ajTj = Kq yi < j < N 

we also have conservation of the quantity 

Nevertheless, the previous quantities may be not well-defined, not even formally. 
Actually, for general configurations of exactly parallel vortex filaments there are 
infinite. In the next section we will bypass this difficulty by the use of renormalized 
quantities. 

4.1.2. Some existence results. First rigorous results on vortex filaments were ob- 
tained by Klein, Majda and Damodaran [39] for pairs of filaments 



+ aiFia^^-i + Fa 
idt^2 + a2T2dl^2 - Fi 



^1 - ^2 

*1 -*2p 



= 



= 0. 



(4.1) 



l^-l - *2|'- 

It may be convenient to introduce the new variables 'J = 4*1 — ^'2, $ = + 4*2 
and Ki = (aiFi + a2F2)/2, K2 = (oiFi — a2F2)/2. Then 



idt^ + K^dl^ + K2dl-^ + (F2 - Ti)^ = 



idt^ + Kidl^ + K2dl<^ + (Fl + F2) 



0. 



(4.2) 
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In [33] (see also [SH1E7]) only the case of identical vortex core parameters ai = 
is considered. Then for identical circulations Fi = r2 we get K2 = so (4.2 1 can be 
decoupled into two equations for <I> and . The solutions have the general wave-like 
form 

where $ is any solution to the linear Schrodinger equation, and where the param- 
eters A, k and w satisfy a consistency equality. On the other hand, for opposite 
circulations we have ki = 0. By Section [3) the exactly parallel solution is given 
by a uniformly translating vortex pair. For nearly parallel filaments one looks for 
solutions with the same symmetry features, i.e. \l/2(i, cr) = — 4'i(t,cr) for all (t, cr). 
The equations (4.2) then yield a coupled system of equations for the real variables 
5Re(*i) = -^12 and 5to(*i) = -i$/2. 

From a mathematical point of view, one can also consider the case when aiFi = 
a2r2 



kq. It follows that K2 = so the equations (|4.2|) are decoupled 



= 



0. 



(4.3) 



In 



For identical circulations one retrieves the previous co-rotating filament pair 
the case of opposite circulations, ^ satisfies the linear Schrodinger equation. 

The linear stability of system (4.1), with ai — around the exactly parallel 
solution to the point vortex system ( 1.4 ) is analyzed in |39j . Linear instability occurs 
when riF2 < 0, while linear stability holds for rir2 > 0. Moreover, numerical 
computations predict global existence in the first case, and self-similar collapse 
in finite time in the second case. So, in contrast with the 2-dimensional setting 
(pairs of vortex points) , finite time collapses of two nearly parallel filaments might 
happen. Following Zhakarov [66j, Klein, Majda and Damodaran suggested that 
such collapses should be well described by a class of self-similar solutions to the 
system (4.2). Also the equations (4.3) could maybe provide a class of self-similar 
collapses. To our knowledge, such constructions have not been rigorously achieved 
yet. 



The next decisive results regarding ( 1.3 1 were accomplished by Kenig, Ponce and 
Vega [37] in 2003. The article [33 is presented in the survey [S5]. In [37], existence 
of vortex filaments in the form of i/^-perturbations of the exactly parallel solution 
Xj(t) is investigated with the ansatz 



'i!j{t,a)^Xj{t)+Uj{t,a), with e C(M, i7^(M)). 



(4.4) 



Note that, thanks to the Sobolev embedding H'^{n.) C L°°(M), the ansatz ( jXij ) 
guarantees that the filaments are well-separated as long as the quantity 
is sufficiently small with respect to the distance between the point vortices Xj{t). 
As long as this remains true, the Uj solve a system of Schrodinger equations with 
potentials that are Lipschitz functions in the Uj : 



i<i<N. 



\Xj - Xk 



- Uk 



\X,-Xk\' 







(4.5) 

This makes it possible in [37] to obtain local well-posedness for (4.5) for any iV, 
for any set of circulations Tj, and for sufficiently small ||Mj(0)|j//i, and at least up 
to times of order min(T*, |lnX]j ll^*j(0)||_f/i |), where T* is the first colhsion time 
between the point vortices. 



12 



On the other hand, Kenig, Ponce and Vega [32 obtained global existence for 
sufficiently small ||uj(0)||//i in the two following special cases: N — 2 and rir2 > 0, 
or A'' = 3 and {Xi{t), X2{t), X3{t)) is the uniformly rotating equilateral triangle 
solution (see Section [s]). The proof of [37j is based on the use of a suitable notion of 
energy. We recall that for filaments given by the ansatz (4.4 1 the quantities defined 



in ^ 4.1.1 do not make sense. This is the reason why |37j introduces the renormalized 
quantities 



A 



T 



\Y.a,T] J Id^M'.l' da- ^- 
T,{\^,\'-\X,\') da 



3 
3 



da 



^k\~\X, 



x,\') 



da. 



Note that in the framework of (4.4) we have \'9j{t,a)-'^k(t,a)\ -)■ \Xj{t) - Xk{t)\ 
when a — )■ ±oo. Moreover, in view of the properties of the point vortex system, the 
previous quantities are still formally preserved in time provided ctjTj — kq for all 
1 < j < iV. 



Finally, we also introduce the quantity 



1(0 = 



jy^k 



\Xj-Xk\{ty 



-l] da. 



which is not necessarily constant, and we define the energy 

s{{^j{t))j) = n + i{t) 



^/Er.r.(-ln 



\^,{a)-^k{a)\' 
\X,-Xk\^ 



\^3~^k\{tr 

\x,-Xk\{ty 



- 1 



da. 
(4.6) 



From now on we will sometimes write 

when not misleading. 

We insist on the fact that, in general, the energy is not constant in time. 

Now, in the cases under consideration in [37] the distances \Xj{t) ~ Xk{t)\ — d 
are constant in time and all the same. Therefore it turns out that 

£{t) = n + ^T = £{0) 

is constant. On the other hand, for circulations having all the same sign it was 
noticed in [37j that the Sobolev embedding combined with the convexity inequality 
(a; — 1)^/2 < X — 1 — Inx < 10{x — 1)^ for x € [3/4, 5/4] imply coercivity for the 
energy: 



— Xlr 



1 



(4.7) 



as long as the vortex filaments ^'j are not too far from the straight filaments Xj 
Finally, it can be easily seen that 



f((X,+y,),)<C^| 



-'ill HI 



(4.8) 
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Therefore in the cases of small iJ^-perturbations considered in [37] the energy is and 
remains small for all time, hence the vortex filaments remain well-separated for all 
time and global existence follows. We stress that the vortex filaments constructed 
in this way remain uniformly close to the straight filaments for all time. 

4.1.3. A class of weak solutions. To conclude this section, we mention that Lions 
and Majda [47] introduced and proved the global existence of "very weak" solutions 



to (1.31 for any > 1, identical circulations and identical vortex core parameters. 



Such solutions satisfy a weak formulation of ( 1.3 ) with suitable test functions, there 
are L-periodic with respect to a, they belong to C(M, H^{0, L)) and they have finite 

energy: supj>Q J2j^k lo I l^j(^)~^fc(OI d-o^ < Actually, such a definition does 
not exclude the possibility of collisions, since the only information is that 

VieE, «'j (t,CT) ^ *fc(t,cr) for a.e. cr. 

Moreover, uniqueness in this class is not known. 

4.2. Symmetric configurations. We present now the results obtained recently 
in [5j . In all the following we will take identical vortex core parameters 

aj = 1 yi<j<N. 

We start with a natural remark that should in part motivate the kind of perturba- 
tions we shall consider then. As we have seen in the previous sections, the dynamics 



of both systems (1.3) and (1.4) is complicated when N becomes large, and in fact 
additional symmetry conditions are needed. In this paragraph, we shall first de- 
scribe a local in time result that extends the one in [3 7) to energy-type spaces. 
We will turn then to the case of four filaments for which global results can be 
obtained without requiring too many symmetries on the perturbations. Finally, 
for any N > 2 we will present a global existence result for dilation-rotation-type 
perturbations of the regular N-polygon exactly parallel configuration. In this latter 
case the perturbations, at any height cr, have a regular N-polygon shape. Moreover, 
traveling waves will be displayed in this setting. We shall complete this subsection 
with another class of simple perturbations. 
For a positive uj we define the energy 

£if) = l [ \daf\'da+'^ [ (|/p-l-ln|/p) da. (4.9) 



2 . 

The motivation for this definition will become clear in the next paragraphs. Re 



current tools in which are partially inspired by I37j (in particular see (4.7) and 



(4.8)), are the following facts (see Lemmas 2.1-.2.3 in [6 ). 

• If / has small energy £(/), then / is close to 1, 

III/P-i||l~<^- 

• If / is small in norm and if t small enough, 

||e"^^/-/|U^< J. 

• If lll/P — 1||l°° < 1/4 then we can compare the energies: 

scpif) ^ Iwd^fwh + ^iii/P - M\h < m < 5fGp(/). 

• If / has small energy £{f) and h is small in H^, then the energy £{f + h) 
is finite. More precisely we have, for absolute numerical constants C,C', 

E{f + h)< CSgpU + h) < C (1 + E{f)) (1 + \\h\\l^) , 
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and 

|||/ + /^|-1||L~ <1. 

4.2.1. N > 2, a local existence result in the energy space. As recalled in §4.1. 2[ local 
existence was proved in [37] for (1.3) for any > 1 and small Tf^-perturbations 
Uj^o = ^j.o ~ ^j,o- The relevant energy functional in that setting is 

So = sii^M 



which has been defined in (4.6 1. We assume next that all the circulations have 
the same sign, for example Tj > for all I < j < N. Then the energy £{{'iij)j) 
satisfies the properties listed above for each function / = — '^^kl/l^j — Afc|. It 
is then natural to consider the problem of perturbations with small energy as a 
single assumption. On the one hand for small initial energy we obtain that ^'o.j 
is not too far from Xj^ and that So < ll'^i.olllfi- On the other hand for 

_o small in it is easy to see that So < C'X^j IKi.olllfi- Therefore the energy 



smallness assumption is weaker than the smallness assumption in general. For 
example, it allows for a larger geometric class of perturbations, namely the small 
rotations. Indeed, take a rotation and translation type perturbation of the square 
configuration, 

with (fo e and ||r'^||//i = 0{e). It follows that So = O(e^) while J2j l|wj,ol!l^i > 
0(1). The idea to use the energy space rather than is also reinforced by the 
fact that in the setting of (1.3), that we shall see later to be related in some cases 
to the Gross-Pitaevskii equation, the L^-norm of the solutions might grow in time 
while the energy remains controllable. 

We next denote by > the minimal distance between the vortices Xj{t) for 
all time (recall that Tj > OVj). In 6J it has been shown the following local well- 
posedness result: 

For an initial configuration (Mj.o)j with small energy, there exists T > and a 
unique solution u — {uj)j G C{[0,T], {M.))^ to the system (4.5) with 

sup \\uj{t)\\Hi < ||ui,o||//i + ^, I <j < N, 

0<t<T 4 

and for T sufficiently small such that 

T(l + So + Y, IKoWm) > (r,),)- 
j 

The proof is based on finding a find a fixed point in the Banach space 

BT = \w^{wi,...,WN)eC{[0,TlH^)'' , sup \\w{t)\\Hi<^ 

L 0<t<T 4 

for the operator A(w) = {Aj{w))j defined by 

. /■* V- r X, + e'^i^^ Q'uj-.o + w, ~Xk- e^"'^''^'tXfc,o - Wk _ X, - Xk \ 
'"^'^h \\X,+e^^^^slu,,o + w,-Xu-e^^''9lu^^o~Wk\^ \X,-Xu\'j 
Then the solution is given by 

Uj[t) = e'*^J'^'uj,o + Wj{t). 

Notice here that we extend Uj locally from a time to not by a fixed point for 
perturbations of the initial data directly, but by a fixed point argument for small 
perturbations Wj of the linear evolutions of the initial data. We must use 
crucially the fact that the deviation e^'^^^Uj^ — Uj,o can be upper-bounded in in 
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terms of the energy at the initial time Sq. This is insured from the previous study of 
the Gross-Pitaevskii equation (see Lemma 3 in [21]). We recall the short proof: the 

Fourier transform of e''*'^'^Uj,Q — Uj^Q can be written as - — ^^j,oiOj the 
norm is bounded by C\/t\\daUjfi\\]^2 and the norm is bounded by C||(?o-Mj.o||L2, 
i.e. 

Ile^^^^'w,, - Uj,o\\„i < C{1 + Vt)\\d„u,^o\\L^ < C{1 + Vt)£o. 
As a consequence of the local in time result we obtain that the solution {uj)j 



to (4.5) exists as long as the energy £{t) remains small. Indeed note that on the 
one hand the norm \\uj{t)\\Hi can grow exponentially, but it cannot blow up as 
long as the energy is sufhciently small. On the other hand, as long as the energy 
remains small, the filaments 'i!j{t,a) remain close the the straight ones Xj{t), so 



no collapse occurs for (1.3 1 



4.2.2. N — A, local and global results around the square configuration. In this part 
we are dealing with {Xj)j the square configuration with equal circulations Tj = 1. 
Let (uj.o)j G i?^(M)* and set 4'j,o = Xj^q + Uj-.q. We introduce the quantity 

^ 11^^1,0 + "3.011^2 , ||U2.0 + M4,o|li2 

So = max <^ So; — — — H ■ — 



In Theorem 1.2. of it is proved that 

// So is small enough, there exists an absolute constant C > 0, and there exists 
a time T , with 



T>Ci 



1 1 



C 1/4 II II 1/2' p 1/3 

So rna,Xj^k\\U],o - Uk fill L2 So 



such that there exists a unique corresponding solution to (1.3) on [0,r], sat- 

isfying '^j = Xj + Uj , with Uj G C ([0, T], 77i(M)), and such that 



3 \^,{t,a)-^kit,a)\ 5 

4^ \X,it) - Xkit)\ ^4' ^^[0'^]' '^^^^ 
Moreover, if the initial perturbation is parallelogram-shaped, namely 

+ U3,o||l2 = ||W2,0 + U4,o||l2 = 0, 

then the solution is globally defined. 

The proof goes as follows. In view of the local existence argument exposed in 
the previous subsection, the extension of a solution is ensured by the control of the 
energy. In this particular setting of the square, direct computations lead to the 
formula 

s{t) = n + -r-A ' + + + 



2 2 



If (^fi, ^'2, *3, *4) is a solution of (1.3 1 then {-'^3, -^4, -^i, -"^2) is also a so- 
lution. So if the initial perturbation is a parallelogram like the one presented in 
Figure [1] it will remain so at later times, hence 

\\{u,+u,){t)\\l. = \\{u2+U4m\\l.=0 

and global existence follows. On the other hand, for a perturbation without sym- 
metry conditions the growth of -|- M3)(i)|P -I- ||(m2 + U4)(i)|P up to time T is 
controlled by a lengthy computation based crucially on the particular properties of 
the underlying square configuration (Xj)j. 
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Figure 1. = 4, Parallelogram-shaped perturbation around the 
rotating square configuration. 



Sit) = -n+i--A+- (hi wiii. + 11(^.2 + u.mwl.) 



Finally, let us notice that one could try the same approach for other configura- 
tions {Xj)j. Explicit computations give for instance 

3 , 3 
2-^+4 

for = 3 and {Xj)j at the ends and the middle of the segment, and 

2 3 

Sit) = +1 - ^ ^ UK- + u,+2 + u,+i){t)\\l, + + ^J+3)it)\\h 

i=i i=i 
for = 6 and iXj)j a regular hexagon. Nevertheless, these quantities have no 
reason to be conserved, unless the perturbations retain the same shape as iXj)j^ 
which enters precisely the framework of the next paragraph. 

4.2.3. N >2, global results for regular polygon-type configurations. In order to get 
global results for a large number N of filaments, in [5^ symmetry conditions arc im- 
posed both on the straight filaments configuration and on the perturbations. More 
precisely, we let iXj)j be the regular TV-polygon configuration, with or without its 
center, rotating with constant velocity lj and centered at the origin. We consider 
filaments i^ j)j with equal circulations T j = 1 and of the type 

■^,{t,<j)^X,{t)^it,a). (4.10) 

So in particular for any time t and any height a, the points jit,a))j form the 
same polygon as iXjit))j up to a rotation or/and dilation, see Figure [2j 

A straightforward computation shows that if j)j is a solution of ( |1.3| ) then $ 
is a solution of 

z9t$ + a2$ + ^JL(i-l$|2)^o. (4.11) 



Eq. (4.11) is an Hamiltonian equation, with Hamiltonian given by 
j j (|<i>|2_i_ln|$|2) da. 

In our context, as we expect that {'^j)j are small nearly parallel perturbations of 
(Xj)j, we have |$| ~ 1, so that (4.11 ) is formally similar to the already mentioned 



Gross-Pitaevskii equation 

iat$ + 9^$ + tj$(i _ |$|2) = 0, (4.12) 
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Figure 2. iV = 6, Rotation and dilation of the polygon rotating 
configuration without center: cr) = Xj{t)^{t, a) for all j. 



with corresponding Hamiltonian given by 

1 
2 



(l-fl 



1 



This formal comparison turns out to be valid: as we have reminded at the beginning 
of this subsection, when |$| ~ 1 the energies f ($) and £gp{^) are comparable. So 
our objective is now to solve globally the Cauchy problem for the equation (4.11) 
in the energy space, in the same spirit as previous works on (4.121, see 

We next assume that w > 0, which is the case as soon as the circulation Tq of 
the central vortex satisfies TV — 1 + 2ro > (see Section [s]). Theorem 1.1 of [6] 
asserts that 



For an initial data 



with {Xj)j a regular polygon configuration, with or without its center, rotating with 
constant positive speed uj, and f ($o) small, (1.3 I has a unique global solution of 
type ( |4.10[ ) such that £{<^{t)) = £($0) and $ - $0 e C(M,i7Hl^))- 1^ particular, if 

$o(cr) ' 1 then for all time '^j{t,a) ' - — >■ Xj{t), so the filaments remain 
nearly parallel. 

The proof is based on the transposition of the arguments presented in §4.2. 1| to 
the setting of (4.11). We notice that 

Since £($) is conserved it follows that the corresponding local existence result for 
(4.11) can be iterated and global existence is achieved. 



In the third part of [6] subsonic traveling waves for (1.31 are constructed, still 
under the condition a; > 0. More precisely, for c < \/2uj sufficiently close to \/2uj 
there exist solutions of (1.3) of the shape 

*j(t,cr) = Xj{t)v{a + ct), 

with v{cr) a smooth function of small energy, of even increasing modulus, exponen- 
tially increasing to 1 as ct tends to 00. Moreover, 

>±C30 



v{a) 



J0± 



< Cv/2w-c2. 
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These traveling waves are the equivalents objects of the grey solitons for the 1-D 
Gross-Pitaevskii equation. The proof of the existence follows the approach in [55] 
and [57] ■ The equation for the profile rj ~ 1 — \v\'^ 

/ \ 1/2 

7] ^ ( - (c^ - 4cj)?7^ + 4:Uj{{r] - 1) ln(l -rj) -7])] 

is a little more complicated than the one for the Gross-Pitaevskii equation which 
can be solved explicitly. It is integrated in 6 using the smallness assumption on 
c — y/2uj- 

As long as the energy is small, the modulus of ^{t,a) stays close to 1 and the 



nonlinearity in (4.11 ) enters the frameworks of [35] (see also [I^), so these traveling 
waves are orbitally stable. Also, the results in [SJj are valid and insure that there 
are no non-trivial supersonic traveling waves. 

Finally, by reusing an argument from [37) , we can transform these traveling waves 



in the following way. System (1.3) is Galilean invariant, which means that for any 
real and any solution the family *j,y(t,cr) = e^^*""'^'''"^ ^j{t,a ~ 2tv) is 

also a solution. By choosing u = ^Juj we get as another solution of ( |1.3[ ) 

(t, a) = e'V^-+^T y(^a + t(c - 2^)), 

which represents a stationary (6*+ — 6'_)-twisted iV- helix filament configuration with 
some localized perturbation traveling in time on each of its filaments. 

4.2.4. TV > 2, shifted perturbations. As noticed in 6 , if the perturbations are 
considered of the shifted form 

^,{t,<j)=Xj{t)+u{t,<j), 

for any {Xj)j with Tj all the same, it follows that u must solve the linear Schrodiger 
equation. For an initial data in L^, by using the dispersion inequality for the linear 
Schrodinger equation it follows that the filaments remains separate globaly in time 
and the perturbations spread along the straight configuration {Xj)j. By considering 
less regular initial data one can get examples of perturbations u decaying at infinity 
that lead to a L°° dispersive blow-up for the linear Schrodinger equation. For 
instance the homogeneous data \x\~p with < p < 1 yields a self-similar linear 
Schrodinger solution, smooth for positive times ([H])- So this solution leads to 



solutions of (1.3) blowing-up in L°° in finite time at height a = 0. Recently an 
example of initial data in but not in that provides dispersive blow-up was 
given in [T^]. More precisely, the solution with initial data e*'^' /(I -I- jxp)™ with 
1/2 < m < 1 blows-up in finite time at one point. This makes a second example of 



solutions of (1.3) that blow up in finite time at a certain height but for which no 
collision occur. 

4.3. Collapses of filaments. In the last part of this paper we describe configura- 
tions of nearly parallel filaments evolving towards collision in finite time. Most of 



them are new results. They are based on perturbations of type (4.10) of exactly par- 
allel filament configurations {Xj)j. Therefore the collapse for the filaments (^j)j 



solutions of ( 1.3 ) is linked to solutions of (4.11 ), with modulus initially close to one 



that vanish at least at one point in finite time. So a pointwise control of solutions 



of (4.11 ) is needed, which is quite unusual in the study of the Schrodinger equation. 

Let us first notice that if (5'j,o)j" leads to a collision in finite time, the shifted 
perturbations introduced in j ]4.2.4[ 

*j-o((7) ^ 1'j-o + uo(ct), 

with uq in for instance, yield a collapse of the same kind. So collapses are not 
isolated phenomena. 



19 





Figure 3. N — 2, Gaussian perturbation around the fixed equi- 
librium solution {Xi,Xq ~ 0,X2 = — ^i), leading to collision in 
finite time. 



4.3.1. Gaussian collapses about a stationary regular polygon configuration. We start 
by recalling the single collapse situation described in [5]. We let {Xj)j be the sta- 
tionary regular A'^-polygon configuration, with circulation Tj = 1 in the N vertices 
of the polygon and circulation —{N — l)/2 in its center (see Sj3|. Note that this 
implies N > 2. Since the angular velocity is a; 0, it follows that $ is a solution 
of the linear Schrodinger equation. The evolution of a real Gaussian e""^ by the 
linear Schrodinger equation can be computed explicitly 

e i+4,t 



VTT4E' 



Therefore we are insured that the non-vanishing initial condition 



^/T^i 



yields a solution for the system (1.3), 



^jit,a)^X,{t) 1 



e i-4i(i-t) 



with '^j—Xj E C (M, (M)) , such that all filament collide together at time t = 1 at 
(7 = 0. Figure |3] displays such initial perturbations for N = 2 around the stationary 
configuration {Xi,X2 — — -^i) with center at Xq = 0. 

As a new observation, we can exhibit a similar scenario of collision for special 
perturbations of the previous data, namely for all initial data of type 



vE',,oM=X,(0) 1- 



with uq small enough in and such that 

e''^uo{y)dy = 
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This is due to the fact that 



g l-4i(l-t) 

v/l-4z(l-t) 



uo(ct) 



is then a solution of the system ( 1.3 1. In view of the dispersion inequahty and the 
explicit formula for a linear Schrodinger solution, we have 

C 

< — 



Vt + i 



■"oIIli, 



uo{0) = C / e' s uo{y)dy. 



Therefore we are insured that at initial time the filaments are separated and 

they collide at time t = 1 at cr = 0. 

4.3.2. Self-similar collapses about a stationary regular polygon configuration. Again 
we consider {Xj)j to be the stationary regular polygon configuration, with circu- 
lations Tj = 1 in the vertices of the polygon and circulation —{N — l)/2 in its 
center. This time we take advantage of the explicit self-similar solution of the lin- 
ear Schrodinger equation exhibited by Cazenave and Weissler [T^ . The initial data 
considered there is ipi^) = l^^l"^ with < p < 1. By combining (3.2), Proposition 
3.3. and Corollary 3.4. b) of [l9] it follows that 



r(V) 



In particular wc obtain the existence of two times < i2 < such that 



< 



1 



2' 

Therefore we consider as initial data the non-vanishing function 



o(a) = X,{0) (l - zSe'*i92V((T)) 



It yields as a solution for (1.3) 



■fj{t,cr)=Xj(Q) (l -iie*(*i-*)aS^((j)) , 



which vanishes at time t 



ti — t2 aX (J — 0. Moreover, Proposition 3.7. in [T^ 

|cr|— foo 



insures that e ""ifi belongs to some space, so we have indeed ^j{t,a) 
X,{t). 

4.3.3. Self-similar collapse around a three-vortex collapse. The purpose of this para- 
graph is to investigate the behavior of vortex filaments around a self-similar col- 
lapsing solution to the three- vortex problem. We consider the one constructed in, 
e.g., [T] or [S3]. The initial configuration is 



-1, Xi(0) = -1, X2(0) = l, ^3(0) = l + i\/2. 



The sufficient collapse conditions P1P2 -I- P2P3 + PsPi = and T — Q are both 
satisfied. Then if we denote by c = — 1(1 + the corresponding center of 

vorticity, the configuration satisfies 



Xj{t)^i{Xj{t)~c)uj{t), l<j<3. 



where 



u:{t) = 



1- i 

T 

The solution is explicitly given by 



Lj = a + ib = 



r = — = 



2b 



3 

71' 



1 



t 



exp — iar In 1 



te[0,r). (4.13) 
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Figure 4. N — 3, Shifted perturbation around a self-similar col- 
lapsing triangle configuration: 'i>j o{a) — Xj{0) + uo{a) for all j. 



As already mentioned, since the exact parallel solution {Xi{t) , X2{t) , Xs{t)) col- 
lapses in finite time then any shifted perturbations of the form 'iij,o{<j) = Xj(Q) + 
Mo(cr), for all 1 < j < 3, with uq € H^, leads to collision at time r (see Figured). 

We next look for other kinds of less trivial perturbations around the triangle 
collapse. As we did before in Section \A.2\ we look for a filament solution under the 
form 

^j{t,a)-c = {Xj{t) - c)<^{t,a), with |$(0,cr)| ~ 1. 
From now on wc will assume that 



Then the equation for $ is 
so finally 



kq, 



1 < i < 3. 



-iXj $ 
X,-cW 



l$l 



0, 



^^t<^> + Ko^l<^> + u;{t)—{l-\'^>\^)^o 



(4.14) 



In the case of an equilateral uniformly rotating triangle we would have io{t) — lo (z 
real and constant and we would retrieve Eq. (4.11). 



We next investigate the existence of self-similar solutions to (4. 14 1. We seek for 
a solution in the form 



(4.15) 



$(^,a)=r(^)exp(^/3(^))exp(^^-^), r(t) > 0, -f{t) eR, 
where r(0) = ro ~ 1 satisfies rg < 1. 



Plugging (4.15) into (4.14) we find the system (recall that rb — 1/2) 

1 



< 



By setting 



2Knr 1 , 

'^+ — + - V\ 

7 r 2(t — t) 

7 - 4ko = 

-r/3+-(l-r2)— =0. 

T T — t 



7(i) = 4Ko(i - ar), 







(4.16) 
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where a > is a parameter, we can find an explicit formula for r: 



a 1 



t 

ar 



t e /, 



^ r ' ^ CXT ' 7 

where / is the largest interval such that / C [0, min(T, ar)) and r > on /. 



(4.17) 



We rewrite (4.17) as 



1 



1 



(4.18) 



and observe that the discriminant in the numerator is negative for a. < 2rQ and 
positive for a > 2r^. We are led to the study of different cases which exhibit 
different kinds of motion. From now on C will denote a positive constant depending 
only on r and a, which can possibly change from one line to another. 

• We have < a < 2r^. It follows that / = [0, min(T, ar)) and 

1 C 

< r{t) < ^^^=^^^= t e [0, min(r, ar)). 



1-^) (1- 



After finding /3 in ( |4.16 ), and taking into account ( |4.13 1 we obtain that the filament 
solution constructed this way satisfies: 

^=<|(*,-c)(t,a)|< ^ 



C 



1- ^ 



Cjl- 



< 



1 - ^ 

ar 

c 



a e 



1 



t 

ar 



t e [0, min(r, ar)). 



If a < 1, this means that a blow-up of the filament solution at each level cr, without 
collapse, takes place at time t = ar (strictly before the occurrence of the three- 
vortex collapse if a < 1 and exactly at the same time if a = 1). If a > 1, no 
blow-up nor collapse occurs for the filaments at the time t = t oi the three-vortex 
collapse. Actually the modulus of the filaments remain regular up to t = r and 
only the angular velocity becomes singular at time r. 



We have 2rn < a < t^tt^-^. 

U — 2(l-r^) 



1 



We formulate ( |4.18[ ) as 
1 



with <, 



2ar2 (1- ^) 
and t 



{t-u){t-t*), tel, 



ar 



1 



Since both t^, > t and 



> T we get the same behavior as for the case 1 < a < 2r§ for the filaments. 

• We have ^{i-r'^) ^ Then > t but < < r; therefore, / = [0,**) is 
strictly contained in [0,t). We conclude that the filament solution evolves toward 
the following self-similar collision in finite time: it exists up to time t* < r, collapses 
at each level a at time t — t* and satisfies 



\^,{t,a)~c\^C\l- 



t 



C\ l- 



t 



t* " ' " V t* 

In particular the collision between the filaments takes place before the collision 
between the point vortices. 



Remark 1. When tq > 1 the last argument does not apply and we are not able to 
construct a collapse in that way. 



23 



Remark 2. In the case of a uniformly rotating polygon (see [4-^^ we have ui{t) = u) 



is constant in time, real, so that the analog of (4.16) is 




Thus the corresponding filament solution is given by 



Uq ( , (lot t 2 ^ \\\ ( 



= X,(t)^exp y — iX - -Y^.riX - -)j j cxp ^^^^^^ _ 

for some uq G S""^ and r G M. It blows up at time t = t but does not collide. 
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